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THE  MICROWAVE  APPLICATIONS  THEORY  PROGRAM  AT  NRL 
AND  SOME  CHEMISTRY  CODE  APPLICATIONS  TO 
IONOSPHERIC  HEATING  BY  MICROWAVE  RADIATION 

1.  INTRODUCTION 

The  advent  of  high  power  pulsed  microwave  devices,  the  magnetrons,  at 
NRL,'*'  which  currently  generate  ^  1  G  Watt  at  X  -  10  cm  and  a  pulse  duration 
of  't  30  nsec,  provide  a  new  opportunity  for  research  in  various  areas  of 
microwave  applications. 

Two  areas  of  interest  in  the  high  power  microwave  applications  are: 
the  microwave  interaction  with  air  (air  breakdown)  and  the  microwave  inter¬ 
action  with  matter  (conductors  and  dielectrics) .  These  two  areas  have  po- 

2 

tential  applications  of  considerable  value,  however,  they  do  require,  in 
general,  a  large  amount  of  power.  The  advent  of  the  magnetron  and  the  free 
electron  lasers  and  their  future  developments  could  provide  the  necessary 
power  using  fewer  devices  compared  to  the  old  ones  with  power  levels  of  10^ 
watts  or  less.  However,  in  order  to  delineate  and  advance  various  applications 
of  high  power  microwave  devices,  a  concerted  effort  in  the  areas  of  air 
breakdown  and  interaction  with  matter  is  desired.  Especially  since  none  of 
the  above  applications  have  been  carried  out  for  microwave  powers  above  10^ 
watts.  To  this  end  the  high  power  microwave  applications  research  at  NRL  is 
conducted  at  two  fronts;  experimentally  and  theoretically.  The  theoretical 
aspect  of  the  program  and  its  goals  are  described  in  this  report.  Furthermore, 
the  report  presents  calculations  for  threshold  power  of  air  breakdown  and  its 
dependence  on  the  microwave  frequency,  pulse  duration  and  the  ambient  pressure. 
In  addition,  some  calculations  are  given  for  the  heating  of  the  ionosphere  and 
the  subsequent  emission  of  the  heated  region. 

Manuscript  submitted  June  16,  1980. 
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2.  THE  MICROWAVE  THEORY  PROGRAM  AND  ITS  GOALS 


« 


The  theoretical  program  can  best  be  described  as  a  basic  research  with 
applications  overtone.  It  is  divided  into  two  categories,  the  atmospheric 
breakdown  and  the  interaction  with  matter. 

2.1  THE  ATMOSPHERIC  BREAKDOWN 

To  transmit  and/or  to  focus  pulsed  microwave  radiation  from  the  ground 
to  any  point  in  space,  the  radiation  must  propagate  through  the  atmosphere. 

The  attenuation  of  the  microwave  radiation,  as  it  propagates  through  air, 
depends  on  various  parameters,  among  these  are:  the  atmospheric  conditions 
(air  density,  temperatures),  the  radiation  frequency,  the  power  density,  the 
pulse  width,  etc.  At  power  levels  near  threshold  for  air  breakdown,  the 
radiation  will  be  blocked  from  reaching  its  destination  due  to  air  breakdown. 
Therefore,  it  is  essential  that  we  understand  fully  the  breakdown  phenomena 
and  its  dependence  on  all  those  parameters.  However,  the  understanding  of 
the  breakdown  to  generate  a  disturbed  atmospheric  condition  may  facilitate  the 
utility  of  air  breakdown  for  certain  useful  applications.  A  disturbed  air 
generated  by  microwave  breakdown  will  emit  radiation  from  uv  to  ir  which  may 
have  some  impact  on  sensors  and  detectors.  The  generated  plasma,  depending 
on  its  size  and  duration,  may  also  affect  communications. 

Therefore,  the  goal  of  the  theoretical  program,  in  this  area,  is  to 
develop  a  theoretical  model  (a  computer  code)  which  describes  the  phenomena 
of  the  radiation  absorption,  the  subsequent  avalanche  ionization,  the  heating 
and  dissociation  of  air  molecules,  the  generation  of  shocks,  overpressures, 
the  emission  of  the  disturbed  air  and  its  subsequent  relaxation.  The  model 
should  incorporate  hydrodynamics  and  a  detailed  air  chemistry  into  a  code  to 
describe  the  phenomena  satisfactorily  for  all  altitudes  and  microwave  fre¬ 
quencies  . 
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The  generated  code,  will  be  utilized  to  predict 

(a)  the  dependence  of  air  breakdown  on  microwave  power  density, 
the  microwave  frequency,  atmospheric  pressure  (altitude), 
the  radiation  pulse  width  and  the  pulse  shape. 

(b)  the  dependence  of  breakdown  on  two  incident  frequencies,  two 
microwave  or  one  microwave  and  an  optical  frequency.  The 
impact  of  two  stage  ionization  and  multiphoton  ionization. 

(c)  shocks,  overpressures  and  their  dependence  on  the  incident 
microwave  power  density,  ambient  pressure,  radiation  pulse 
shape  and  width. 

(d)  spatial  and  time  dependence  of  the  air  species,  degree  of 
ionization,  dissociation,  electron,  ion  and  cluster  ion 
densities  (depending  on  altitudes,  etc.)  and  various  kinetic 
temperatures  e.g.  the  electron  temperature,  the  vibrational 
temperature  and  the  neutral  temperature. 

(e)  emission  from  the  disturbed  air  (uv,  visible  ir,  etc.).  Its 
magnitude,  duration  and  its  dependence  on  microwave  power 
density,  frequency,  pulse  width,  pulse  shape  and  the  altitude. 

(f)  plasma  absorption  length,  power  per  unit  volume  expended  for 
breakdown  to  generalize  towards  larger  volumes.  Power,  number 
of  pulses  and  their  time  separation  needed  to  sustain  such  a 
plasma. 

(g)  relaxation  of  the  disturbed  air  and  the  impact  of  the  late  time 
air  chemistry  on  multi  pulse  breakdown  characteristics. 

Another  interesting  aspect  is  the  air  breakdown  and  its  characteristics 
near  surfaces.  This  phenomenon  will  also  be  studied  and  provides  a  link  to 
the  interaction  of  radiation  with  matter. 

2.2  MICROWAVE  INTERACTION  WITH  MATTER 


This  aspect  of  the  theoretical  program  concentrates  on  the  interaction 
of  microwave  radiation  with  conductors  and  dielectrics  in  vacuum  and  rarefied 
atmosphere.  Its  purpose  is  to  model  the  physical  changes  in  matter  and  how 
they  are  brought  about  by  the  incident  radiation?  These  physical  changes 
are  the  heating  of  the  target  (its  temperature),  melting  and  vaporization. 

And  how  these  changes  depend  on  incident  power  densities,  radiation  wave¬ 
lengths,  the  pulse  width  and  the  pulse  shape?  The  theoretical  model  which 
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describes  this  phenomena  entails  radiation  absorption,  heat  equation  vapor¬ 
ization,  breakdown  in  vaporized  region,  the  hydrodynamics  of  the  vaporized 
region,  energy  coupling  and  decoupling  to  the  solid  matter  by  the  vaporized 
and  ionized  region.  These  basic  elements  are  to  be  incorporated  into  a  com¬ 
puter  code  which  will  be  utilized  to  predict: 

(a)  the  temperature  of  the  matter  under  irradiation,  the  melting 
point  and  the  vaporization.  The  dependence  of  these  changes 
on  the  microwave  power  density,  microwave  frequency,  the  dura¬ 
tion  of  the  pulse  and  the  shape  of  the  pulse  and  the  angle  of 
incidence. 

(b)  breakdown  of  the  vaporized  atomic  species,  its  hydrodynamic 
development,  degree  of  ionization  and  the  range  of  emission. 

(c)  the  effect  of  the  broken  down  vapor  on  coupling  or  decoupling 
the  microwave  energy  to  the  solid  matter. 

(d)  the  emitted  radiation  from  the  vaporized  and  ionized  region  and 
their  absorption  into  the  matter. 

(e)  generated  shocks  or  detonation  waves  and  their  impact  on  the 
matter. 

It  is  obvious  from  these  brief  statements,  on  the  goals  of  the  theoret¬ 
ical  program, • that  large  computer  codes  with  detailed  ab  initio  type  calcula¬ 
tions  are  needed.  However,  estimates  can  be  made,  using  simple  theoretical 
models,  to  gain  an  understanding  of  some  or  more  aspects  of  the  problem. 

These  approaches  will  be  utilized  as  one  proceeds  to  develop  the  detailed 
computer  model  which  in  principle  should  provide  completeness  to  the  solution 
of  the  problem. 

Simple  theoretical  models  are  pursued  in  the  next  section  to  obtain  some 
characteristics  of  the  air  breakdown  phenomena. 
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3.  AIR  BREAKDOWN 

In  this  section  a  brief  theory  of  the  microwave  air  breakdown  is  given. 

Quantitative  relations  are  derived  for  the  breakdown  threshold  power,  its 
dependence  on  the  microwave  wavelength  and  pulse  lengths.  Calculations  are 
presented  for  various  altitudes. 

3.1  THEORY  OF  BREAKDOWN 

The  breakdown  in  air  is  exemplified  by  the  generation  of  a  plasma.  This 

occurs  when  sufficient  power  from  a  microwave  source  is  focussed  in  air, 

whereupon  the  free  electrons,  present,  gain  energy  under  the  influence  of  the 

intense  electric  field  associated  with  the  focussed  radiation.  The  energy 

gained  by  electrons  is  expended  in  producing  more  electrons  by  electron  impact 

ionization,  culminating  in  cascade  ionization. 

There  have  been  considerable  experimental  studies  of  breakdown  and 

ionization  rates  in  many  gaseous  elements  and  in  air  under  the  influence 

3-10 

of  microwave  radiation.  With  the  advent  of  lasers  and  high  power 

Q-switched  lasers,  however,  the  air  and  gaseous  breakdown  studies  at  optical 

11  12 

frequencies  have  increased  considerably.  ’ 

The  mechanism  for  gas  breakdown  can  be  explained  by  the  classical 

4 

microwave  breakdown  theory.  The  electron  gains  energy  from  the  electric 
field  only  when  it  collides  with  an  atom  and  the  sudden  change  in  its  velo¬ 
city  results  in  the  transfer  of  the  oscillation  energy  gained  from  the  field 

13 

into  random  motion.  This  energy  gain  process  has  been  shown  to  correspond 
to  the  quantum  mechanical  description  of  the  energy  gain  by  the  free  elec¬ 
trons  from  the  radiation  field  through  the  free-free  transitions  in  the  field 
of  an  atom. 

The  equation  of  motion  of  a  free  electron,  in  air,  under  the  influence 
of  an  oscillating  electric  field  is 

I 
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m  —  =  -  v  mv-eE  exp(-iwt)  (1) 

at  m 

where  is  the  collision  frequency  for  momentum  transfer,  oj  is  the  frequency 

of  the  radiation  field,  E,  and  v  is  the  velocity  of  the  electron.  In  Equa- 

-► 

tion  (1),  the  term  v  mv  represents  a  continous  resistive  damping  force  on 

m 

the  electron  due  to  collisions  with  atoms  and  molecules.  The  real  part  of  the 
velocity,  in  phase  with  the  oscillating  field  is 

ev  Exp(-iu)t) 

_  ni  / \ 


2  2 
m(v  +  oj  ) 
m 


The  rate  of  energy  gain  by  the  free  electron 


^  J  v.E.  Exp(-iu)t)dt  ^  is 


de  \  _  e2E2  _ ji 

dt  1  mv  2 

/  g  m  uj  + 


2  2 
2  E  v  4 

where  E  *  — r - is  the  effective  electric  field  (equivalent  to  a  static 

e  00-  +  v  2 

m 

field) . 

On  the  other  hand,  the  electron  loses  its  energy  through  elastic  and 
inelastic  collisions.  Electron  energy  is  also  lost  by  attachment,  recombina¬ 
tion  and  diffusion.  The  total  energy  lost  can  be  expressed  as 


-  \  D  £  2m  / 

—  -  (Vi 


+  v£  +  Zve+v 

a  e  e  e  t 


where  D  is  the  diffusion  coefficient,  A  the  diffusion  length,  e  is  the  average 
electron  energy,  v . ,  V  ,  v  and  v  are  the  collision  frequencies  for  ioniza- 
tion,  attachment,  excitations  and  effective  recombinations,  respectively,  while 
and  eg  are  the  corresponding  energies  for  ionization  and  excitations. 
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respectively.  The  second  term  on  the  right  hand  side  of  Equation  (4)  repre¬ 
sents  the  energy  loss  by  electrons  through  elastic  collisions.  Therefore,  the 
net  energy  gained  by  electrons  is 


d e 
dt 


(5) 


Furthermore,  the  time  development  of  the  electron  density,  n^,  can  be  des¬ 
cribed  by 


dn 

e  D 

-t—  =  nv  -nv  -n  —  -  n  v 
dt  ei  ea  e  .2  er 


(6) 


where  it  is  obvious  that  electron  attachment  and,  diffusion  and  recombination 

will  control  the  cascade  ionization,  and  must  be  overcome  for  breakdown  to 

occur.  Two  limiting  cases  of  Equation  (6)  can  be  of  great  interest.  These 

are  the  diffusion  controlled  and  attachment  controlled  cases  of  air  breakdown. 

At  atmospheric  pressure,  however,  the  air  breakdown  is  attachment  controlled. 

To  estimate  the  magnitude  of  the  terms  in  (6),  consider  the  condition  for  a 

/dne  \ 

cw  breakdown  I  -j—  =  01  which  implies  that 


v 


i 


v 

a 


(7) 


The  diffusion  coefficient  can  be  expressed  as 

D  p  -  3.2  x  105E  (8) 

where  p  is  the  pressure  in  Torrs  and  C  is  the  average  electron  energy  in  units 

of  eV.  Thus  for  typical  energies  of  few  eV  and  pressure  of  one  atmosphere 

2 

(760  mm),  the  diffusion  coefficient  is  D  -  8  *  10^  *  This  implies  that 

2  3-1 

the  rate  of  diffusion  from  an  area  of  1  cm  is  10  sec  which  is  very  small 
compared  to  the  attachment  rate  of  electrons  to  0^.  There  are  two  distinct 
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attachment  processes  that  an  electron  undergoes  in  air.  These  are:  the 
three-body  attachment  to  02  i.e. 

e  +  02  +  02  -*•  0^  +  02  (9) 

e  +  02  +  N2  -*■  0~  +  N2  (10) 

and  the  dissociative  attachment 


e  +  02  -*■  0  +0  (11) 

14 

The  rate  coefficients  for  reactions  (9)  and  (10)  are 

3.5  *  10  Exp(-0.052/T  )  and  10  cra^/sec,  respectively.  At  room 

e 

temperature,  the  attachment  rate  due  to  reactions  (9)  and  (10)  is 

^  6  x  10^  sec  1  which  is  larger  than  the  diffusion  rate  by  many  orders  of 

magnitude.  Reaction  (11),  however,  is  highly  dependent  on  the  electron  energy 

where  the  cross  section15  is  appreciable  for  electron  energies  (4.5-10  eV) , 

-2  2 

the  peak  of  this  cross  section  is  10  tt  aQ  at  electron  energy  of  6.5  eV. 
Recombination  in  air  proceeds  through  dissociative  recombinations 

N2+  +  e  -*■  N  +  N  (12) 

02+  +  e  -  o  +  0  (13) 

and  by  two  and  three  body  recombinations.  These  can  be  presented  by  the 
following  relations: 

A+  +  e  -►  A  +  hv  (14) 

A++e  +  e-*-A  +  e  (15) 

A++e  +  M-»-A  +  M  (16) 

•’here  A+  is  an  atomic  or  a  molecular  ion  and  M  is  a  third  body.  These 
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recombinations  (Reactions  12-16)  are  all  electron  temperature  dependent  and 
their  coefficients  are  given  in  Table  1,  taken  from  Reference  16. 


TABLE  1 

Recombination  Rate  Coefficients 

T  in  (eV) 
e 


Reaction 

Rate  Coefficient 

12 

,  -  -8,  .-0.39 

3 

-1 

4.3  x  10  (T  ) 
e 

cm 

sec 

13 

_8  -0.63 

3 

-1 

2.1  x  io  ° (T  ) 
e 

cm 

sec 

14 

-12  ° 

3 

-1 

^  10  at  300  K 

cm 

sec 

15 

-27  -4  5 

5.2  x  10  L  (T  ) 
e 

6 

cm 

-1 

sec 

16 

-31  -2  5 

6.5  x  io  Ji(T  ) 
e 

6 

cm 

-1 

sec 

c 

The  critical  electron  density,  Ng,  for  a  plasma  generated  by  radiation 
with  a  wavelength  A  is 

c  13-2 

N  -  10  A  (17) 

e 

which  implies  that  for  microwave  radiation  with  A  >  1  cm  the  critical  electron 
c  13  “8 

density  is  Ng  4  10  cm  For  such  a  condition  and  for  a  pressure  of  one 
atmosphere  the  recombination  is,  generally,  dissociative  (see  Reactions  12 

and  13)  especially  when  electron  temperatures  are  'v  few  eV.  Thus  for  Tg  -  1  eV, 

13  —3  —ft  13 

N  ■  10  cm  the  dissociative  recombination  rate  is  4.30  *  10  x  10 
e 

^  4.3  x  10^  sec  ^  which  is  15  times  smaller  than  the  attachment  rate. 
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3.2  BREAKDOWN  THRESHOLD  POWER 

The  air  breakdown  calculation,  in  principle,  can  be  carried  out  by 
solving  Equation  (5)  and  (6)  with  the  knowledge  of  the  electron  velocity 
distribution  and  the  appropriate  atomic  processes  whose  rates  are  electron 
energy  dependent.  However,  reasonable  estimates  (lower  limits)  can  be  ob¬ 
tained  for  breakdown  threshold  power  using  approximate  solutions.  The  rate 
of  energy  growth  for  the  electron  in  the  radiation  field.  Equation  (3),  can 
be  written  in  terms  of  the  incident  radiation  intensity,  I,  as 


de 

dt 


2 

4rre  I 
me 


v 


m 


0.11  I  v 


m 


2  2 
0)  +  v 

m 


2  2 

Od  +  V 

m 


(18) 


E  X  2 

where  —  =  —  was  utilized  (E  =  27  /i  (V/cm)  with  I  in  watts/cm  .  The  thres¬ 
hold  condition  is  obtained  using  Equation  (18)  assuming  that  energy  gained 
by  the  electron  is  expended  to  ionizations  and  that  other  loss  mechanisms  are 
ignored.  Thus  the  energy  gained  by  the  electron  during  the  radiation  pulse. 
At,  is  equal  to  the  ionization  potential  of  the  atom  or  the  molecule,  1^, 
times  the  critical  number  of  generations,  Kcr, 


At 

f  —■  dt  -  K  I. 

dt  cr  i 

o 


Thus  the  threshold  power  for  air  breakdown  is 


(19) 


th 


2  2 
I .  K  (u)  +  v  ) 
i  cr  _ m 

0.11  At  v 

m 


(20) 


Some  interesting  observations  can  be  made  of  the  result  in  Equation  20. 
First,  the  threshold  power  for  breakdown  increases  inversely  with  the 
radiation  pulse  width  i.e.  if  one  could  generate  a  very  short  pulse  then 
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one  may  transmit  higher  powers  in  air  without  breakdown.  Second,  for  the 

2  2 

microwave  frequency  go  «  v  ,  the  threshold  power  increases  as  V  i.e.  with 

m  m 

pressure  and  is  independent  of  the  microwave  frequency  (high  pressure  regime) . 

For  a  pulsed  breakdown  the  cascade  electron  density  build  up  can  be 
described  by 


N£  -  Nfi(0)  eVAt  =  Ne(0)  2  cr  (21) 

Where  v  is  the  net  ionization  rate  and  Ng(0)  is  the  initial  density  of  the 

_3 

free  electrons,  generally  ^  1  cm  .  Furthermore,  it  is  assumed  that  at  the 

end  of  the  radiation  pulse  the  electron  density  reaches  its  critical  value, 
cr 

(see  Equation  17),  where 


K 


cr 


(22) 


To  estimate  the  threshold  power  for  air  breakdown  under  optical  fre¬ 
quencies,  Equation  (20)  can  be  utilized  with  Kcr  values  ranging1"*"  from  36  to 
40.  Using  K  ■  40  and  1^  =  14.8  eV  for  the  average  ionization  potential  of 
air  molecules  in  Equation  (20)  results  in 


th 


8.6  x  10 


-16 


u>2  + 


m 


At  v 

ra 


(watts/cm  ) 


(23) 
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The  collision  frequency  can  be  estimated  to  be  v  =  5.3  x  10  p  where  p  is 
the  air  pressure  in  Torrs,  therefore,  Equation  (23)  reduces  to 


8.6  x  io-16  co2  +  (5.3  x  109  p)2 

xth - At - ; :  ~  9 — 

5.3  x  10  p 


(24) 


For  optical  region  go  =  101A  which  is  much  larger  than  v  ,  one  obtains  for  air 

ID 

10  2 

at  atmospheric  pressure  a  breakdown  threshold  power  of  2  x  10  (watts/cm  ) 
using  a  pulse  with  duration  of  10  ^  sec.  This  value  is  of  the  same  order  but 


11 


I 


« 


» 
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lower  by  a  factor  of  three  compared  to  the  breakdown  threshold  power  obtained 
18 

experimentally,  indicating  that  the  classical  approach  can  explain  the 
breakdown  process.  Thus  the  breakdown  calculations  require  the  interrelation 
of  many  factors  which  are  the  ionization  rate,  the  energy  loss  rate,  the 
average  electron  energy,  the  electric  field  and  finally,  the  threshold  power 
for  breakdown.  This  requires  a  more  detailed  analysis  which  is  carried  out 
in  the  next  section  for  microwave  breakdown  of  air. 


3.3  AN  ANALYTIC  METHOD  FOR  THRESHOLD  POWER  CALCULATION 

Various  approaches  to  air  breakdown  calculations  exist. ^  Some 

start  with  the  knowledge  of  the  ionization  rate  of  air,  V.,  which  generally 
Vi 

is  expressed  as  —  .  This  quantity  depends  on  the  electron  energy  which  in 
P  E 

turn  depends  on  the  critical  electric  field,  —  .  For  a  Maxwellian  electron 

P 

velocity  distribution  and  for  an  average  electron  energy  of  e  <  20  eV,  the 

Vi  20 

following  relation  for  —  has  been  obtained 

P 


—  =  2.5  x  io 
P 


7  [’• 


6  (e)*5  +  0.6  (e)3/2 


]  M 


(25) 


where  it  is  assumed  that  in  air  ionization  is  basically  from  oxygen  molecule 

20 

(I,  -  12.06  eV) .  Furthermore,  an  empirical  relation  has  also  been  found 
E 

—  Q 

between  e  and  —  where 
P 


e  -  0.0675  — 
P 


(26) 


This  relation  utilized  into  (25)  yields 


which  is  shown  in  Figure  1.  Thus  the  problem  is  a  simple  one  for  air  break¬ 
down  when  it  is  attachment  controlled.  For  cw  the  condition  for  breakdown  is 


12 


(28) 


or  equivalently 


V 

a 


(29) 


v  E 

ci  6 

The  value  of  —  is  known  and  thus  one  can  obtain  —  from  Equation  (27) .  How- 
P  P 

ever,  for  pulsed  breakdown,  the  ionization  rate  must  be  high  enough  so  that 
the  electron  density  reaches  its  critical  value  within  the  duration  of  the 
radiation  pulse.  This  requires  that 


v 

a 

P 


e 


N  (0) 
e 


(30) 


or 


v.  v 

i  a  1  . 

p  p  pAt  °®e 


I"l013/A2  •  Ne(0)l 

L  V.  J 


(31) 


Thus  for  any  wavelength.  At  and  p  we  can  obtain  — .  From  Figure  1  or  Equation 

E  P 

£ 

(27)  we  obtain  the  corresponding  —  and  hence  the  threshold  breakdown  power. 

Using  A  *  3  cm,  p  *  760  mm.  At  ■  10  ^  and  N  (0)  -  1  in  Equation  (31)  we  obtain 
E  6 

—  =  30  volts/cm  Torr. 

P 


An  application  of  Equation  (31)  with  Ng(0)  ■  1  for  a  fixed  pulse  duration 
(At  *  10  ^sec)  is  shown  in  Figure  2  for  several  different  wavelengths.  In 
this  Figure  — ' s  are  given  as  a  function  of  air  pressure  (altitude).  The 
corresponding  effective  electric  fields  as  a  function  of  pressure  (altitude) 
are  shown  in  Figure  3.  In  this  figure  the  independence  of  the  electric  field 
on  the  radiation  wavelength  with  increasing  pressure,  is  a  manifest  of  the 
attachment  controlled  domain  for  breakdown. 
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3.4  AIR  BREAKDOWN  CALCULATION  WITH  DIFFUSION 

The  calculation  in  Section  3.2  was  carried  out  ignoring  the  diffusion 
effect.  However,  when  the  density  of  air  is  reduced,  diffusion  may  become  an 
important  factor  in  air  breakdown  and,  therefore,  it  has  to  be  considered. 

In  this  section  we  present  air  breakdown  calculations  by  including  the 
diffusion  effects  phenomenologically.  In  addition,  we  present  the  effect  of 
the  radiation  pulse  width  on  the  breakdown  thresholds. 

The  criterion  for  a  cw  breakdown  is 


For  a  pulsed  condition,  however,  the  criterion  for  breakdown  is 


(32) 


F  fi  -  *a  -  %)  ■  m  l°*e  (C'V0')  <33> 

This  implies  that  two  other  parameters  have  to  be  known  in  order  to  calculate 
the  critical  field.  These  are,  D  and  A.  The  diffusion  coefficient  is  related 
to  Ee/P  according  to 

DP*  ^29.0  +  0.9  -jfj  x  104  (34) 

and  A  is  generally  a  characteristic  length  of  the  breakdown  geometry  and  can 
be  estimated.  In  Figure  4  we  present  breakdown  threshold  power  calculations 
as  a  function  of  pressure  for  several  wavelengths  of  the  incident  radiation. 
These  calculations  were  carried  out  with  A  ■  0.75  cm  and  At  *  3  *  10  ^  sec. 

It  is  obvious  from  this  Figure  that  the  threshold  breakdown  power  increases 
with  increased  frequency  of  the  incident  radiation.  However,  these  calcula¬ 
tions  were  repeated  for  a  shorter  pulse  width  i.e.  At  *  3  x  10  7  sec.  to 
demonstrate  the  effect  of  shorter  pulse  duration.  These  latter  results  are 
shown  in  Figure  5  where  it  is  obvious  that  the  reduction  in  the  radiation 
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pulse  width,  for  each  wavelength,  results  in  higher  breakdown  threshold 
power. 
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4.  SOME  APPLICATIONS  USING  AN  NRL  CHEMISTRY  CODE 

21  22 

The  NRL  Master  Code,  ’  an  E  and  F  regions  chemistry  code,  was  utilized 

for  the  ionospheric  electron  heating  and  increased  ionization.  The  code  is 

21  22 

described  in  detail  elsewhere.  *  However,  a  brief  description  of  the 
code  is  in  order. 

The  code  solves  a  set  of  coupled  rate  equations  which  predicts  the  time 
histories  of  the  following  species: 

N2,  N2(A),  N2+,  02,  02(a1A),  02(b1i:),  02+,  02+(a47T),  NO,  N0+,  N,  N(2D), 

N(2p),  N+,  N+(1D),  nVs),  0,  0(1D) ,  0(XS),  0+,  0+(2D),  0+(2P),  ne,  Tg,  Ty 
and  T  ,  where  for  the  species  the  ground  state  is  implied  whenever  no  elec¬ 
tronic  state  is  indicated.  The  electron  temperature  T  ,  the  vibrational 

temperature  of  the  nitrogen  molecule,  T^,  and  the  kinetic  temperature  of  the 

20  21 

heavy  particles,  T  are  also  calculated.  *  The  code  can  accept  any  ioni- 

3 

zation  and  heating  source  to  disturb  the  ionosphere  and  it  is  purely  a  chem¬ 
istry  code  and  does  not  contain  transport  or  hydrodynamic  effects,  however, 
it  can  be  coupled  with  such  a  code  for  a  complete  description  of  the  dis¬ 
turbed  air  phenomena.  The  code,  furthermore,  can  calculate  the  emission  from 
various  excited  states  of  the  ionospheric  species.  It  can  be  used  to  gain 
further  insight  into  microwave  deposition  into  the  ionosphere. 

We  have  selected  an  altitude  of  250  Km  to  deposit  microwave  radiation 

_2 

with  A  -  3  cm  and  a  pulse  duration  of  At  ■  10  sec.  Three  sets  of  incident 

3  2 

power  densities  were  used  1,  10  and  10  Watt/cm  and  some  selected  results  of 

these  calculations  are  shown  in  Figures  (6-8).  At  a  power  density  of 
2 

1  Watt/cm  no  change  occurs  in  the  ambient  electron  density  (Fig.  6).  The 
energy  gained  by  the  ambient  electrons,  however,  is  expended  in  inelastic 
collisions  resulting  in  excitations  of  various  electronic  states.  The  time 
histories  of  two  such  states,  0(^S)  and  0(^D)  of  the  oxygen  atom  are  shown 
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23 

in  Figures  7  and  8.  They  are  metastable  states  with  life  times  of 

^  0.6  sec  and  135  sec,  respectively.  They  emit  radiation  at  X  -  5577  A  and 

o  24 

the  6300  A,  respectively,  and  are  the  well  known  auroral  green  and  red 

oxygen  emissions. 

3  2 

When  the  incident  power  density  is  raised  to  10  and  10  Watts/cm  , 

however,  the  electron  heating  and  multiplication  occurs  resulting  in  peak 

5  -3  5  —7 

electron  densities  of  ^  2  x  10  cm  and  8  x  10  cm  ,  respectively.  For 

these  cases  the  enhanced  0(^S)  and  0(^D)  densities  are  shown  in  Figures  7  and 

2  3  5  -3 

8  with  peak  densities  at  10  Watts/cm  being  8.5  *  10  and  2  *  10  cm  , 

3  7  5  -3 

respectively.  The  respective  values  at  10  Watts/cm  are  1.5  x  10  cm  and 
1.7  x  io7  cm-3. 

4.1  THE  MICROWAVE  ABSORPTION  COEFFICIENT 

In  the  previous  section  some  results  of  the  energy  deposition  in  the 
ionosphere  were  presented.  It  is  of  interest,  however,  to  know  the  penetra¬ 
tion  depth  or  the  attenuation  coefficient  for  the  microwave  radiation  in  air 
or  gaseous  medium  in  general.  This  quantity  is  essential  for  various  appli¬ 
cations  e.g.,  the  propagation  length  and  the  size  of  a  plasma  which  can  be 
generated  by  a  radiation  of  a  given  intensity  and  wavelength. 

The  process  of  radiation  absorption  by  the  free  electrons  is  the  inverse 
Bremsstrahlung  with  two  components.  These  are  due  to  the  electron  neutral 

and  the  electron  ion  collisions  in  the  field  of  an  electromagnetic  radiation. 

18  25 

The  corresponding  absorption  coefficients  are  ’ 


a  -  10”3A  n  N  \2  ST 
n  e 


(35) 


and 


ot 

P 


2  x  10 


-23 


,2 

n  N  a 


e 


(36) 
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respectively.  In  Equations  (35)  and  (36)  T  is  the  electron  temperature  in 

degrees  Kelvin,  N  and  N+  are  the  neutral  and  the  ion  densities,  respectively 

and  are  valid  for  regions  of  hv  «  kT.  Furthermore,  in  Equation  (35)  the 

-15  2 

momentum  transfer  cross  section  is  assumed  to  be  ^10  cm  .  It  is  clear 
from  these  relations  that  when  the  degree  of  ionization  is  low,  the  dominant 
absorption  mechanism  is  the  inverse  neutral  Bremsstrahlung  and  when  the  de¬ 
gree  of  ionization  rises  the  electron  ion  component  becomes  important  and  in 
many  cases  will  predominate,  especially  when  the  breakdown  is  due  to  optical 
frequencies  i.e.  (lasers). 

From  these  relations,  it  is  obvious,  that  for  our  applications  of 
microwave  radiation  at  250  Km  with  electron  density  of  CLO3  -  107)  cm  3  and 
neutral  density  ^  109  cm-3,  the  absorption  length  is  very  long  as  if  the 
radiation  will  never  be  attenuated.  Therefore,  one  can  generate  a  column  of 
excited  states  many,  many  kilometers  in  length  which  will  emit  the  green  and 
the  red  oxygen  lines.  As  an  example  the  peak  emissions  of  S  and  D  for  an 
incident  microwave  power  of  1000  Watt /cm2  are  2.5  *  105  and  1.2  *  105  photons 
per  cm3,  respectively.  For  a  column  of  1  Km  the  emission  from  the  green  and 
the  red  lines  will  be  300  kR  and  144  kR,  respectively. 

4.2  IONOSPHERIC  MODIFICATION  AND  ITS  IMPACT  ON  COMMUNICATION 

In  the  previous  section  we  gave  estimates  of  the  microwave  power  needed 
to  increase  the  ionization  level  in  the  ionosphere.  An  obvious  impact  of 
Che  increased  electron  density  would  be  the  disruption  of  communications,  at 
higher  frequencies,  between  satellites  and  the  earth.  The  disruption  occurs 
when  the  ionosphere  acts  as  a  mirror  for  the  transmitted  frequency.  This 
frequency  is  related  to  the  critical  electron  density  by  Equation  (17). 

More  detailed  studies  on  this  area  will  be  forthcoming. 
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5.  AIR  BREAKDOWN  SPECTRA 

19  -3 

The  air  at  the  sea  level  is  a  mixture  of  ^(2.1  x  10  cm  )  and 
18  —3 

02(5.38  *  10  cm  )  i.e.  with  a  ratio  of  4  to  1.  The  molecular  densities 

26 

fall  off  exponentially  with  increasing  height.  There  occurs,  however,  a 
break  from  the  exponential  fall  off  for  O2  around  100  km,  where  the  dissoci¬ 
ation  of  0^  becomes  important,  giving  rise  to  the  oxygen  atom  which  becomes 

one  of  the  dominant  species  at  higher  altitudes.  For  example,  at  250  km  the 

8  -3  3  -3  26 

N2  and  0  densities  are  2  x  10  cm  and  7.6  x  10  cm  ,  respectively. 

The  air  breakdown  spectra  will  be  due  to  band  and  line  emissions  from 
N2»  N2+*  °2’  °2+’  °*  0+’  N  an<*  N+  an<*  var^ous  other  species  which  are  gen¬ 
erated  as  a  result  of  the  breakdown.  To  this  one  must  add  plasma  emissions 
due  to  various  processes  such  as  free-free  free-bound,  -bound-bound,  etc.  A 
detailed  discussion  on  air  breakdown  emissions  will  be  the  subject  of  another 
report. 

For  air  breakdown  with  microwave  radiation  at  the  sea  level,  however, 
under  a  single  pulse  (10  ^  -  10  ^sec)  operation,  one  may  observe  the  emission 

spectra  of  Nj  and  especially  from  the  second  positive  bands  of  ^  and 

+  27 

the  first  negative  bands  of  ^  .  These  two  band  systems  correspond  to 

+  2  2  3  3 

N2  (B  I  +  X  Z)  and  ^(C  -  B  tt  )  transitions,  respectively.  The  emission 

o 

in  the  second  positive  bands  system  extends  from  2800  -  5000  A  with  its  (0,0) 

o 

transition,  at  3371  A,  being  the  strongest.  The  emission  in  the  first  nega- 

o 

tive  bands  system,  on  the  other  hand,  extend  from  3300  -  5900  A  with  its  (0,0) 

o 

transition,  at  3914  A,  being  the  strongest.  In  Table  II  some  strong  bands  of 

28 

these  systems  are  presented  along  with  their  radiative  transition  rates. 
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TABLE  II 


Strong  Optical  Transitions  From  ^  Second  Positive 


and 

First  Negative  Bands 

Systems 

Transition 

CM) 

0 

Wavelength  (A) 

Transition  Rate 

(0,0) 

3371 

1.11(7)* 

(0,1) 

3576 

7.27(6) 

(0,2) 

3805 

2.83(6) 

(1,0) 

3159 

1.02(7) 

(1,2) 

3536 

4.7(6) 

(1,3) 

3755 

3.93(6) 

(1,4) 

3998 

1.24(7) 

KJ 

4 

o 

o 

3914 

1.24(7) 

(0,1) 

4278 

2.2(6) 

(1,2) 

4236 

2.66(6) 

(1,3) 

4651 

1.93(6) 

(*)  1.11(7)  read  1.11  x  107 
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Fig.  1  —  The  ionization  frequency  per  Torr  v\!p  as  a  function  of  Ee/p 
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2  —  The  ionization  frequency  per  Torr  as  a  function  of  pressure  for 
several  incident  radiation  wavelength  with  At  =  10“ ®  sec 
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Fig.  3  —  The  corresponding  effective  electric  field  for  Fig.  2 
as  a  function  of  pressure 
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Fig.  5  —  Same  calculation  as  in  Fig.  4  for  a  shorter  radiation 
pulse  width  (3  x  10— ?  sec) 
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Fig-  7  —  The  time  history  of  0(^8),  which  emits  the  green  line  at 
5677  A,  at  an  altitude  of  250  Km  for  microwave  radiation  power 
densities  of  1,  10  and  1000  Watt/cm^.  Hie  radiation  pulse  width 
is  IQ-2  sec. 
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0('D)  DENSITY  (cm-3) 


Fig.  8  —  The  time  history  of  0(^D),  which  emits  the  oxygen  red  line, 
at  6300  A,  at  an  altitude  of  250  Km  for  microwave  radiation  power 
densities  of  1,  10  and  1000  Watts/cm^.  The  radiation  pulse  width 
is  10“2sec. 
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